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Abstract. We consider the conditions under which the Cayley trans- 
form of the Kronecker product of two Hermitian matrices can be again 
presented as a Kronecker product of two matrices and, if so, if it is a 
product of the Cayley transforms of the two Hermitian matrices. 
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1. Introduction 

Let Mn be the algebra of all n x n matrices over the complex field and 
Hn C Mn the subalgebra of Hermitian matrices. As usual, a conjugate 
transpose of a complex matrix A £ Mn will be denoted by A* . Now, suppose 
that A G Hn, i.e., A* = A, and let In be the n x n identity matrix. Then 
{A + iln)~^ exists and 

UA = {A-iIn){A + iIn)-^ 

is called a Cayley transform of A. It is easy to see that Ua is a unitary 
matrix and the inverse transform is given by 

A = i{In + UA){In-UA)-y 

Furthermore, +1 cannot be an eigenvalue of Ua- In the following we give 
some basic examples. 

Example 1.1. 

(1) If ^ = /„, then Ua = -Hn- 

(2) If A is the n x n zero matrix, i.e., A = On, then Ua = —In- 

(3) If ^ is a diagonal matrix, then Ua is also a diagonal matrix. 

(4) If A has degenerate eigenvalues, then Ua has degenerate eigenvalues 
as well. 

(5) If y4 G Hn is unitary, i.e., = In, then Ua = —iA. 

The Cayley transform is actually a generalization of a mapping of the 
complex plane to itself, given by 

U{z) = ^^, zeC\{-i}. 
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In particular, U maps the upper half plane of C conformally onto the unit 
disc of C and the real line ffi injectively into the unit circle. Moreover, no 
finite point on the real line can be mapped to +1 on the unit circle. 

Let us continue with some useful properties of the Cayley transform. 

(1) liV €Mn is invertible, then UyAV-^ = VVaV'^ for A G i/„. 

(2) If X S C" is an eigenvector for an eigenvalue A G M of a matrix 
A G Hn, then x is an eigenvector of a Cayley transform Ua and 
U{\) = (A — i)/(A + i) is its eigenvalue. 

(3) li Ae Hm and B G H^, then Ub(^a = PUa®bP\ where P G Mmn 
is the permutation matrix satisfying P{A® B)P^ = B®A. Here, P* 
denotes the transpose of a matrix P and (g) denotes the Kronecker 
product (see, for example, [9l flO]). 

(4) liA,Be Hra such that [A, B] = AB - BA = then [Ua, Ub] = 0. 

The Cayley transform is named after Arthur Cayley (see (21 E])- In the 
last few decades, a lot of results about the Cayley transform and its ap- 
plications, mostly in mathematics and physics, have been obtained. For 
example, Calixto and Perez-Romero [1] utilized the Cayley transform for a 
complex Minkowski space. Jadczyk [5] applied the Cayley transform in the 
compactification of the Minkowski space. Furthermore, Eisner and Zwart 
[1] studied Co-semigroups and the Cayley transform. Noncommutative Cay- 
ley transforms have been introduced by Popescu [7] and an application of 
the Cayley transform for rotation of elasticity tensors has been studied by 
Norris [6]. 

In mathematical physics, the main applications of the Cayley transform 
is to the Hermitian matrix 

' a + d b — ic^ 
b + ic a — d , 



H 



with a, 6, c, d real. Here a question is: What is the condition on a, b, c, d such 
that the matrix H and the Cayley transform Uh coincide (perhaps up to a 
phase)? Note that the eigenvalues of H are 



Ai,2 = a± \/W+(F+(f 
and the eigenvalues of Uh are 

\io-i Afo-l „. Aij 



Now, it is easy to see that the eigenvalues of H and the eigenvalues of the 
Cayley transform Uh never coincide since Ai^2 are real (if Ai^2 = 0, then 
H = O2 and Uh = — -^2)- On the other hand, if the eigenvalues of H and the 
eigenvalues of Uh differ only by a phase, then 

A1.2 - ^ _ I X I 
Al,2 + « 
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This yields that Ai^2 = il- In each case the phase difference is —i since 
f^(Ai,2) = —i^i,2, i-e., Uh = —iH. In particular, one of the following holds: 

(1) a = 1 and b = c = d = 0, 

(2) a = —1 and b = c = d = 0, 

(3) a = and a'^ + b^ + c'^ = 1. 

In the paper, we discuss the following question. Let A £ Hm and B G Hn 
be two Hermitian matrices. Then the Cayley transform provides the unitary 
matrices Ua and Ub, respectively. Now, A<^$B is again a Hermitian matrix 
and the Cayley transform gives us another unitary matrix Ua®b- 

(1) What is the condition on A and B such that Ua<sb can he again 
presented as a Kronecker product of two matrices? 

(2) What is the condition on A and B such that Ua®b can he presented 
as a Kronecker product of the Cayley transforms of two Hermitian 
matrices? 

(3) What is the condition on A and B such that Ua®b = Ua ® Ub ? 

To conclude our introduction, we give some examples of Hermitian matri- 
ces for each of the questions above. First we consider two Hermitian matrices 
A and B such that the Cayley transform Ua(^b cannot be presented as a 
Kronecker product of two complex matrices. 

Example 1.2. Let A = B = diag(l, 0) = ( J q ) be diagonal 2x2 Hermitian 
matrices. Then Ua®b = diag(— i, — 1, — 1, — 1) is a diagonal 4x4 matrix 
which cannot be presented as a Kronecker product of two 2x2 complex 
matrices. 

Now, let us write one simple example showing that Ua®b = Ua ®Ub 
does not hold in general. In this example the Cayley transform Ua^b can 
be presented as a Kronecker product of two complex matrices. 

Example 1.3. Let A = O^n be the m x m zero matrix and i? = 0^ be the 
nxn zero matrix. Then Ua = -Im, Ub = -Im, Ua^b = -Imn = Uj^^Uj^, 
and Ua®Ub = Imn- Obviously, Ua(^b ^Ua®Ub- 

Finally, we give four simple examples of pairs of 2 x 2 Hermitian matrices 
{A, B) satisfying Ua(^b = Ua®Ub- 

Example 1.4. 

(1) A = -h and B = {1± V2)l2 

(2) A = -h and B = diag(l ± a/2, 1 =F \/2) 

2. Cayley transform on Hermitian matrices 

Let A G Hm and B € Hn be two Hermitian matrices. In this section, 
we first answer the question (Theorem 12. 2^ : Under what conditions can we 
write the Cayley transform of A® B as a Kronecker product C ® D with 
C G Mm, D G Mn? The next natural question here is (Theorem 12. 4p : If 
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Ua®b = C ®D for some C G Mm, D G Mn, under what conditions C = Ua 
and D = Ub? 

In the following we use a variation of the result about the separability 
constraints which was proved in [8]. 

Lemma 2.1. Let { Ci, . . . , Cm } be a basis for an m- dimensional subspace 
of Mm and { Di, . . . , } be a basis for an n- dimensional subspace of Mn- 
Then a matrix A in the form 

m n 

-4 = ^ ^ Oj- fcCj (g) Dk 
can be written as Ai ® A2 with A\ G Mm and A2 G Mn if and only if 

for all p,r £ {1, . . . , m} and q, s £ {1, . . . , n} . 

Now we are in the position to write our first result. 

Theorem 2.2. Let A G Hm CL^d B G Hn- Then the Cayley transform of 
A®B can be written asC®D G Mm®Mn if and only if one of the following 
conditions is fulfilled. 

(a) A has one eigenvalue. 

(b) B has one eigenvalue. 

(c) A has two eigenvalues 0.1,02 and B has two eigenvalues 61, &2 such 
that aibia2b2 = 1. 

Proof. Let A G Hm be a Hermitian matrix with eigenvalues {ai, . . . ,am} 
and corresponding orthonormal eigenvectors {xi, . . . ,Xm} and let B G Hn 
be a Hermitian matrix with eigenvalues {61, . . . , bn} and corresponding or- 
thonormal eigenvectors {yi, . . . ,yn}. Then the Cayley transform of ^ ® i? 
can be written in the form 

m n , 

UAm = 2^2^ ^T— ^^j^i ® ykVk ■ 
j=i k=i ~^ 

We already know that the eigenvectors are preserved under the Cayley trans- 
form. Thus, by Lemma l2.ll (with identifying Cj — )• xjXj, ykVki ^'^d 
0''j,k — ^ {o,jbk — i)/{ajbk + i)), we find that 

(1) Ua^b = C®D 

for some C G Mm and D G M„ if and only if 

Opbq — i Orbs — i Opbs — i Orbq — i 
Opbq + i Orbs + i Opbs + i Orbq + i 

for all 1 < j,p < m and 1 < k,q < n. This equation can be rewritten as 

(2) (op - ar){bq - bs){aparbqbs - 1) = 0. 
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It follows that either ajbkapbq = 1 or aj = ap or bk = bq must hold for all 
eigenvalues aj, Up and bj,bq of A and B, respectively. 

Obviously, if either (a), (b), or (c) holds, then ([2]) is fulfilled and we are 
done. For the converse implication, suppose that equation ([2]) holds and 
assume that neither A nor B has only one eigenvalue. Let ap,ar and bq,bs 
be two distinct eigenvalues of A and B respectively. It follows from ([2]) that 
apttrbqbs = 1. Thus, Op 7^ 0, 7^ 0, bq 7^ 0, and bs / 0. Let at be an 
eigenvalue of A and assume at is distinct from ap and a^. It follows that 
Qpatbqbs = apQrbqbs = 1 which yields at = ar, a contradiction. Thus, A has 
2 eigenvalues and they are nonsingular. Similarly B has 2 eigenvalues and 
they are nonsingular. If we denote the eigenvalues of A with ai, 02 and the 
eigenvalues of B with 61,62, then the equation ([2]) reduces to 01026162 = 1. 
The proof is completed. □ 

Theorem 2.3. Let A £ Hm and B G Hn be such that the Cayley transform 
Ua»b = C (Si D' £ Mm (S) Mn- Then there exist C G Hm and D G such 
that UA(g)B = Uc®Ud- 

Proof. We may assume that C' and D' are unitary. Furthermore, we may 
assume that neither C' nor D' have 1 as an eigenvalue. Namely, there exists 
G M such that neither e*^C" nor e~^^D' have 1 as an eigenvalue (since C 
an D' are finite). Thus, by the invertibility of the Cayley transform, there 
exist C G Hm and D G Hn such that Uc = C and Ud = D' ■ Consequently, 
Ua®b = Uc®Ud- □ 

The next theorem answers the question when Ua^b = Ua ®Ub- 

Theorem 2.4. Let A G Hm and B G Hn- Then Ua^b = Ua ® Ub if and 
only if one of the following conditions is fulfilled. 

(a) A has one eigenvalue o 7^ and B has one or two eigenvalues given 
by 

, a(l-a)±y^a2(l-a)2-4o 

61 2 = . 

2a 

(b) B has one eigenvalue 6 7^ and A has one or two eigenvalues given 
by 

6(1 - 6) ± V 62(1 -6)2 -46 



ai,2 



26 



Proof. Let A G Hm be a Hermitian matrix with eigenvalues {oi, . . . ,Om,} 
and corresponding orthonormal eigenvectors {xi, . . . ,Xm} and let B G Hn 
be a Hermitian matrix with eigenvalues {61, . . . , 6^} and corresponding or- 
thonormal eigenvectors {t/i, . . . ,yn}- Then the Cayley transform of A, B 
and A® B can be written as 
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III lb 7 

j=l k=l 

Comparing Ua ®Ub and Ua^b yields 

ttj — i bk — i CLjbk — i 



CLj + i bk + i ajhk + i 

or equivalently 

(3) ajbk{l- aj - bk) = 1, 

where 1 < j < m and 1 < k < n. So, Ua^b = Ua (X" Ub if and only if the 
relation ([3]) holds for all 1 < j < m and 1 < A; < n. 

Now, if either (a) or (b) holds, then it is easy to see that ([3]) is fulfilled 
and we are done. For the converse implication, suppose that the relation 
([3]) holds for all 1 < j < m and 1 < k < n. Since we assumed that 
Ua(^b = Ua®Ub, Theorem 12.21 implies three cases. 

Case 1. Suppose that A has one eigenvalue a. Then 

abl - a(l - a)bk + 1 = 

for all 1 < /c < n. Note that a 7^ 0. Thus, it is easy to compute that 



a{l-a)±Ja^il-aY-Aa 

= 2-a ■ 

In particular, B has one or two eigenvalues as in the case (a). 

Case 2. If B has one eigenvalue, then, using the same arguments as in the 
previous case, we obtain (b). 

Case 3. Finally, suppose that A has two eigenvalues ai,a2 and B has two 
eigenvalues 61, 62 such that ai 61 02^2 = 1 and which satisfy Q for all j, k G 
{1,2}. However, no such 01,02,61,62 £ ^ exist. Namely, from ^ we have 
the following four equations 

ai6i(l - (ai + hi)) = 1, 0162(1 - (fli + 62)) = 1, 

0261(1 - (02 + 61)) = 1, 0262(1 - (02 + 62)) = 1. 

Using 01610262 = 1, we obtain 

(1 - (oi + 61)) = 0262, (1 - (oi + 62)) = O261, 

(1 - (02 + 61)) = 0162, (1 - (02 + 62)) = 0161. 

Thus, we find 

1 = 02620161 = (1 - (01 + 6i))(l - (02 + 62)) 
= 02610162 = (1 - (ai + 62))(1 - (02 + 61)). 

It follows that 

(oi - 02) (61 - 62) = 
which cannot be satisfied since oi 7^ 02 and 61/62. Hence, this case does 
not yield any solutions. □ 
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Corollary 2.5. For all A E Hm, we have Ua^a 7^ Ua C?) Ua- 

Proof. If Ua ® Ua = Uai^iA then, by Theorem 12.41 ^ = o/m for some a G 
M \ {0} with 

^ q(1-q) _ _ 4„ = 0, 

2a V ; > 

which has no solutions for a G M. □ 

Corollary 2.6. For all B G Hn, we have Uj^^b 7^ Ui^ (X" Ub- Similarly, 
for all A G Hm, we have Ua^i„ ^^Ua® Ui^. 

Remark 2.7. The natural question here is whether the analogue results 
hold true if we replace the Kronecker product with some other product, for 
example, with the star product or with the direct sum. In particular, if 
A G Hm and B G Hn, then one can easily show that Uai^b = Ua ® Ub- 
Now, recall that the star product of ^ = (ajfc) £ M2 and B G M„ is defined 

by 

(an Oi 
xn 0,12 \ 
Onxl B 0„xl , 
xn O22 I 

where 0,^x1 is a column of n zeros, and Oixn = O^xi- ^'^ ^^^^ known that 
there exists a permutation matrix P G M„,4.2 such that P{A-kB)P* = AqB. 
Therefore, if A G H2 and B G Hn, then UA-kB = Ua *Ub as well. 

Remark 2.8. At the end, let us point out that we have considered just 
the bipartite case, i.e.. Mm ® M„ with integers m,n > 2. But we can 
naturally extend our results to the multipartite systems (8) • • • (8) 
with 77,1, ... , Um > 2 and m > 2 when 

f/M„^®-®M„„ = Um^^ (g) f/M„^^^®-®M„„ Vj G {1, . . . , m - 1} 

or 

UM„^^-(giM„. = Um„^ ® UM^^^i-iSiMn. Vj G {2, . . . , m}. 
Namely, if, for example, m = 3 and Ai G -ffm, A2 G Hn2, A3 G Hn^, then 
we first observe when ?74i®(yl2(g>A3) = Ua^ C/A2(g)A3- According to Theorem 
12.41 this is true if and only if one of the following conditions is fulfilled. 

(a) Ai has one eigenvalue (which is nonsingular) and ^42 (8) A^ has one or 
two eigenvalues (which are nonsingular) given by the exact formula 
(see case (a) in Theorem 12. 4p . 

(b) A2 (8> ^3 has one eigenvalue (which is nonsingular) and Ai has one or 
two eigenvalues (which are nonsingular) given by the exact formula 
(see case (b) in Theorem 12. 4p . 

Again, using Theorem 12.41 we find out that Uai»A2^A3 = Uai ® Ua2 ® Ua^ 
if for distinct j,k,l G {1, 2, 3}, Aj and A^ have one eigenvalue (nonsingular) 
and Ai has one or two eigenvalues (nonsingular eigenvalues are given by the 
exact formulas). Similarly, UAi^-^Am = Uai • • • Ua^ for Hermitian 
matrices Ai G -ffm, • • • , Am G Hn^, m > 3, if all the matrices Ai, . . . , Am, 



8 



YORICK HARDY, AJDA FOSNER, AND WILLI-HANS STEEB 



with one possible exception (this exception may have two eigenvalues, both 
nonsingular), have one eigenvalue (nonsingular) given by the exact formula. 
We omit the details since the proofs are rather technical. 

Remark 2.9. The results in the bibartite case do not extend to the multi- 
partite case in a straightforward way. For example, we have 

Ulm'S)Im'S)Im / ^ Uj^ (g) C//„^ , 

Ulm®Im®Im®Im / ^^/m (g) Uj^ ® C//^ 

but 

m 
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